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Abstract. We compute the shear and bulk viscosities, as well as the thermal 
conductivity of an ultrarelativistic fluid obeying the relativistic Boltzmann equation in 
2 + 1 space-time dimensions. The relativistic Boltzmann equation is taken in the single 
relaxation time approximation, based on two approaches, the first, due to Marie and 
using the Eckart decomposition, and the second, proposed by Anderson and Witting 
and using the Landau-Lifshitz decomposition. In both cases, the local equilibrium is 
given by a Maxwell- Jiittner distribution. It is shown that, apart from slightly different 
numerical prefactors, the two models lead to a different dependence of the transport 
coefficients on the fluid temperature, quadratic and linear, for the case of Marie and 
Anderson- Witting, respectively. However, by modifying the Marie model according 
to the prescriptions given in Ref. pQ, it is found that the temperature dependence 
becomes the same as for the Anderson- Witting model. 
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1. Introduction 

The study of the transport properties of 2D relativistic fluids from the standpoint of 
kinetic theory is an important topic, still awaiting a complete systematization. Kremer 
and Devecchi [2] calculated the bulk viscosity of a two dimensional relativistic gas using 
the Anderson- Witting collision operator and the Chapman-Enskog expansion (3j SJ [5] , 
but did not investigate the shear viscosity and thermal conductivity. 

In this paper, we compute the transport coefficients, namely the bulk and shear 
viscosities and the thermal conductivity, by using two single relaxation time models (also 
called model equations). The first one, proposed by Marie [6], is appropriate for mildly 
relativistic fluids with moderate values of the Lorentz factor, 7 < 2. The second one, by 
Anderson and Witting [3], on the other hand, can deal with significantly larger Lorentz 
factors. The Marie model, as it was initially proposed in Ref. [6], is not appropriate 
to describe a gas of ultrarelativistic particles, due to the fact that it implies an infinite 
relaxation time in the limit where the mass of the particles becomes zero, thus leading 
to divergent transport coefficients. However, Takamoto et al. pQ proposed a modified 
Marie model, whereby the relaxation time of the Boltzmann equation is redefined in 
such a way as to regulate the aforementioned infinities. In addition, it is known that by 
promoting the relaxation time to the status of a dynamic field, it is possible to describe 
complex flows far from equilibrium, such as they occur in turbulence [7j. Therefore, this 
single relaxation time model will also be included in the present study of the transport 
coefficients. In general, model equations do not give the same transport coefficients 
as the ones obtained from the full Boltzmann equation. However, it was proven in the 
Ref. [8] that the methods of Chapman-Enskog and Grad lead to the same approximations 
to transport coefficients, when polynomial expansions of the distribution function in the 
peculiar velocity are performed. 

In both cases, we use the moment expansion of the non-equilibrium distribution, 
similar to the fourteen fields^ HI [5] in the three-dimensional case. To the best of 
our knowledge, this task has never been undertaken before for the case of two spatial 
dimensions. This is all but an academic exercise, since two-dimensional relativistic 
flows arise in many areas of modern physics, say, cosmology, e.g. in galaxy formation 
from fluctuations in the early universe [10], as well as in high-energy nuclear physics, 
e.g energetic heavy ions collisions [llj. Two-dimensional ultrarelativistic fluids received 
a further boost of popularity in 2004, with the discovery of the gapless semiconductor 
graphene [T2l[T3j . This consists of literally a single carbon monolayer and represents the 
first instance of a truly two-dimensional material (the "ultimate flat land" [H]), where 
electrons move like massless chiral particles, whose dynamics is governed by the Dirac 
equation, with the Fermi velocity playing the role of the speed of light in relativity 
[T5l [T6] . However, the calculation of the transport coefficients is more general and 
can be extended to any statistical system of quasi-particles governed by relativistic 
Boltzmann-like equations, i.e. it might apply to a whole class of systems where physical 
signals are forced to move close to the their ultimate limiting speed |17j . 
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The results of this paper are restricted by the range of applicability of 
the Boltzmann equation to (quasi) two-dimensional systems. It is well known 
that linearizing hydrodynamics in two dimensions leads to divergent transport 
coefficients, both in classical and relativistic systems [18j [19]. However, as long as 
the Boltzmann equation provides a useful semi-phenomenological approximation to 
transport phenomena, as for example evidenced by the use of the Boltzmann equation 
in quantum transport [20], results of our computations remain valid. 

We wish to emphasize that the main goal of the present paper is to derive the 
transport coefficients for 2 + 1 dimensional relativistic fluids, out of prescribed relaxation 
times. In the non-relativistic case, this task is pretty straightforward, since in an absolute 
reference frame, there is no ambiguity as to the definition of the macroscopic observables 
(kinetic moments) in terms of the Boltzmann distribution. In the relativistic case, 
on the other hand, this correspondence, i.e. the projection from the kinetic to the 
hydrodynamic space, is much less direct and requires careful consideration. Besides 
its theoretical interest, the practical target of this work is to provide operational input 
for lattice formulations of the Boltzmann equation, which have recently shown major 
potential for the numerical simulation of a broad class of relativistic flows across scales, 
from astrophysical flows, all the way down to quark-gluon plasmas [211 1221 |23j 123], 
including turbulent phenomena in the two-dimensional electronic gas in graphene |25j . 



2. Non-Equilibrium Distribution 

The single relaxation time Boltzmann equation for the Minkowski metric, i] al3 , can be 
written as [I] 

vnc 

P ^f = -—(f-n , (i) 

TM 

for the case of the Marie model [6], and as 



p»U, 



P"dJ = -^(f-n , (2) 

for the case of Anderson and Witting model [3]. Here, m is the mass of the particles, c 
the speed of light, ks the Boltzmann constant, / the probability distribution function 
(which can denote any scalar field in phase space), r the single relaxation time for the 
Anderson- Witting model, and tm the respective one for the case of the Marie model. 
The 3-momentum is denoted by p M = (p°,p), and the macroscopic 3-velocity by 
Greek indices run from to 2, being the temporal component, and we have adopted 
the Einstein notation (repeated indexes are summed). For the purpose of this work, we 
are using the signature (+, — , — ). The equilibrium distribution / cq is given by [I] 

r = A(n,T)eM-p,U»/k B T) , (3) 

where A(n,T) is a normalization constant that depends on the temperature and the 
number of particles density n. In this work, we will study the ultrarelativistic regime, 
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which is characterized by £ = mc 2 jk^T <C 1. From now on, we will use natural units, 
m = c = kg = 1, and the following notation: 

A a/3 = r] af3 -U a U p , 

rp<a(3> _ rp(a(3) _ £ ^ 

V Q = A Q/3 <9 /3 . (4) 

In order to identify the physical meaning of the different terms in the balance and 
transport equations, it is useful to introduce decompositions of these terms with respect 
to orthogonal quantities. Note that A a/3 and U a are orthogonal quantities, A a/3 Up = 0, 
so that any 3-vector can be decomposed into this orthogonal basis. We begin with the 
Eckart decomposition [2B], and later make the due corrections to take into account the 
one proposed by Landau and Lifshitz [31 H] . 

In the Eckart decomposition, the entropy 3-flow, defined by 

d 2 p 

can be written as follows: 



S a = -Jp a fHf)^- , (5) 



S a = nsU a + tp a , (6) 

where s = S a U a /n is the entropy per particle and tp a = A^S 13 the entropy flux. 

In order to obtain the non-equilibrium distribution, we begin by maximizing the 
entropy per particle, under the following constraints: 

,d 2 p 
,d 2 p 



N a U a = U a 



J"? 

«Jp a p?f l - 

T <^>"U a = U a Jp a p^pP>f^ , (7) 



p° 



where, 



and 



p <» p P> = p (« p P) _ -A^^pV' , (8) 



p(V> = i(A^A{ + A?Af)pV . (9) 

In principle, the moments N a and T a ° would be a more natural choice. However, 
the resulting procedure to compute the Lagrange multipliers via entropy maximization, 
while leading to equivalent results [J], proves significantly more complicated. The 
problem of maximizing the entropy is equivalent to consider the following functional, 

jr = s _ \N a U a - \pT a(3 U a - A <7/3> T <7/3>Q ?7 Q , (10) 
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and applying the functional derivative 8T /8f = 0. Here, A, \@, and A <7j g> are Lagrange 
multiplier that we must determine. In two dimensions, there are nine independent 
multipliers, since by definition ?7 7/3 A <7/ 3> = 0, which represents an extra equation. 

Following the procedure, as a result, one can approximate the non-equilibrium 
distribution function by, 

-n(A + A/ + A <7 ^V)] . (11) 
By decomposing the Lagrange multipliers in the orthogonal basis in space-time, 

X^X'Up + X'A} , (12) 



A 



<7/3> = AUpU 7 + ^A a (A^ + A£[/ 7 ) 



+ A 



ad 



-A a5 A 



(13) 



and inserting these variables into the equilibrium distribution, Eq. ( II II) . we obtain 



1 — n\ — n [X'U, 



j, + a;a 7 ) 

n^AU^ + ~A a (AyUp + A%U 1} 



+ A c 



-A tw, A 7/3 ) pV 



(14) 



In the Grad method, we need to determine the value of the Lagrange multipliers 
in terms of the macroscopic fields, n, U a , T, u, p <a ^> ; and q a (being the particle 
density, macroscopic 3-velocity, temperature, dynamic pressure, pressure deviator, and 
heat flux, respectively). The dynamic pressure is defined by u = —fiS/ a U a , and the 
pressure deviator by P <a P> = 2r]'V <a U l3> , being \i and r\ the bulk and shear viscosities 
respectively. In this procedure, we also need the moments of the equilibrium distribution 



function, which have been introduced in Appendix |Appendix A 



Let us impose that the actual distribution function carries the same first moment 
as the equilibrium distribution, namely: 

N ° = S lp "T^l rp "T ' (15) 

and apply the projectors, U a and A^, to obtain respectively the first two equations for 
the Lagrange multipliers, 

- n 2 [X + 2T(3TA + A')] = , 

nTA 7 A 7/3 + 3nT 2 A 7 A 7/3 = . (16) 
In order to obtain the other equations, we calculate the energy-momentum tensor, 

p° 



r* = J /pV^ 



(17) 
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and introduce the projectors: 

A^A* - £A^A a/} ) T Q/3 = P<^> , 

A^T^ = <f , 

A Q/3 T Q/3 = _2( p + w ) , 

f/^T^ =e , (18) 

where e and p are the energy density and hydrostatic pressure, respectively. 

Thus, by inserting the distribution function in Eq. i fPTl) . and taking the projectors 
defined in Eq. ( TT8|) . we obtain the following relations, 

00 = -n 2 T(X + 3T(4TA + A')) , 

e = 2nT — 2n 2 T(\ + 3T(4TA + A')) , 

g f = 3n 2 T 2 A 57 A^ + 12n 2 T 3 A 57 A 5 , 

p<j5> = _ 6n 2 T 3 A<7 5> _ ^ 

Note that by imposing the state equation for the ultrarelativistic system, e = 2nT, 
implies that that u becomes zero. This is equivalent to say that the bulk viscosity 
vanishes, like in three dimensions. In addition, this gives A = — TA', A = — A'/6T, 
where A' can take any value. For simplicity, we will set A' = 0. The arbitrariness on 
this parameter is due to the fact that, in the ultrarelativistic regime (virtually massless 
excitations), the number of particles density, n, and the temperature of the system, are 
not independent, since n ~ T 2 (e.g. gas of photons). In other words, the number of 
particles density is fixed once the energy density e of the system is chosen, and therefore 
one Lagrange multiplier falls apart. 

In order to obtain the other Lagrange multipliers, we solve the system of algebraic 
equations, Eqs. (fTBT) and ( fl9l) . to obtain: 

d 3n 2 T 3 ^ ' 

7 Y< 7 5> = 1 p<y 8 > 
6n 2 T 3 

A^X' s = -^-q^ . (20) 

Replacing these equations into the definition of the non-equilibrium distribution, 
we obtain: 



cq 



(21) 



nT 2 3nT^ ' pP P + QnT 3 P P 
This is the non-equilibrium distribution function for a two-dimensional ultrarela- 
tivistic system, as expressed in terms of the nine moments. The results described here 
can also be obtained by using the so-called triangle scheme [27], which is equivalent 
to the Grad method. In order to calculate the explicit values of the heat flux and the 
pressure deviator, we have to solve the Boltzmann equation. For the purpose of this 
study, we choose two approaches for the collision operator, the first one proposed by 
Marie (6], and the second one by Anderson and Witting [3]. 
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T <*Pi = f fpCpPpyHP. . (22) 

J U 



The third order moment of the distribution function can be calculated as follows: 

By replacing Eq. (l2Tj) into this equation and rising the indexes for the nine fields, 
we obtain for the third order moment, 

nl z 

„er/"A p<e\> 

Note that for an accurate calculation of the third order moment of the distribution 
function, knowledge up to the fifth order moment of the equilibrium distribution 
(denoted by subindex E) is required. The moments of the equilibrium distribution, 



Eq. ([3]), are introduced in Appendix Appendix A Thus, replacing the respective 
moments of the equilibrium distribution, we obtain the third order moment, 

T a ^ = lhnT 2 U a U l3 W 

- 3nT 2 (r] afS U J + rf^U^ + rf^U a ) 

- 2T(r l aP q 1 + rfV + ^V) 

+ 10T([/ a £/V + tWV + U p Wq a ) 

+ 5T(p <a/3> U~ ( + p< a ~i>U p + p <pn<> U a ) . (24) 

With this expression at hand, we are ready to consider the Boltzmann equation. 
For the case of the Marie model, we have all the needed quantities in place. However, 
for the Anderson Witting, some corrections are required, which we shall be introduced 
in Sec. 13.21 



3.1. Marie Model 

For the case of the Boltzmann equation for the Marie model, Eq. ([TJ), it is assumed that, 




These conditions are satisfied in the Eckart decomposition [26J, as in the three 
dimensional case. The first order moment of the distribution, N a , is defined by nU a , 
the second moment is defined by 

and the third order moment is calculated as described before, via Eq. ([21]). These are 
functions of U a , so that we do not need any correction to the moment definitions. By 
integrating the Boltzmann equation, Eq. (pQ), in the momentum space, and taking into 
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account the relations (J25]), we obtain, d a N a = 0, and by multiplying the equation by 
p 13 and repeating the same procedure, we further obtain d a T al3 = 0. These are the 
conservation equations for N a and T al3 . 

However, by multiplying the equations by p^p 7 , we obtain a different equation, 
d a T alBl = — (1/ta/)(T^ 7 — Tj 7 ), which contains the information about the transport 
coefficients. By a standard iteration procedure [I], we can convert this equation into 

Tf*-Tg> = -T M d a T% f * . (27) 

This means that for the Marie model, we just need to know the third order moment of the 
equilibrium distribution, and the second order moment of both, the non-equilibrium and 
equilibrium distributions. The corresponding transport coefficients will be calculated in 
Sec.H 



3.2. Anderson- Witting Model 

For the case of the Anderson- Witting model, we should use the Landau-Lifshitz 
decomposition [3l H]. Such decomposition implies that U a must be calculated, by 
solving the eigenvalue problem, T^U^a = tUia (subindex L denotes Landau-Lifshitz). 
In general, U a , calculated with the Eckart decomposition using N a = nU a , will differ 
from the one calculated with the energy flux, U". As a consequence, we must find the 
relation between both quantities and the correct expression for the third order kinetic 
moment. 

In the Landau-Lifshitz decomposition, we assume that N a U a = N^U a , U a T al3 = 
U a T^P . Moreover, the first and second order moments of the distribution are defined 
by, 

N a =nU1 + J a , 

T°f> = P< a f»- PL A? + e L UZUl . (28) 

It can be easily shown that the correspondence between U a and U£ for the 
ultrarelativisitic case is U a = — q a /3n, with J a = —q a /3T. The conservation 
equations d a N a = and d a T a/s = can be obtained by multiplying by 1 and p 13 , and 
integrating the Boltzmann equation in the momentum space, respectively. 

By multiplying by p^p 1 and applying the Maxwell iteration procedure as before, 
we obtain: 

(T a/37 -Tf J )U La = -rdalf 1 ■ (29) 

Note that in this case, at variance with the Marie case, the third order moment of 
the non-equilibrium and equilibrium distribution functions is needed. To calculate the 
correct expression for the third order moment, it is sufficient to replace U a by Uf i —q a /2>n 
into Eq. ( 12^1) . retaining up to linear terms in the nine fields. This delivers: 

- 3nT 2 { V a[S Ul + rj aj Ul + n^Ul) 
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- T(rfV + ?f V + W) 

+ 5T(U2Ufr + Ufptf + UlUlq a ) 

+ hT(P <aP> Ul + P<^>u[ + P<^>U1) . (30) 

Everything being in place, we next proceed to calculate the transport coefficients 
for the two-dimensional ultrarelativistic system, using both decompositions. 



4. Transport coefficients 

Since we found that the bulk viscosity vanishes, we focus on the shear viscosity and the 
thermal conductivity. First, let us consider the Marie model and use the expressions 
( fi~8l By applying the projector ApU~, to Eq. ( !27l) . we obtain the heat flux, 

q 5 = 3nTr M (v 5 T - i-vV) , (31) 

and by applying the projector ApA^ — |A e<5 A / 3 7 , we obtain the pressure deviator, 

P«*P> = 6nT 2 r M V <a U p> . (32) 

From this two relations we can conclude that the transport coefficients in the model of 
Marie are given by, 

km = 3C ^ B nr M , rj M = ^nk B Tr M , /j m = , (33) 

being the thermal conductivity, the shear viscosity, and the bulk viscosity, respectively. 
Note that we have restablished the physical units. 

For the case of the Anderson- Witting model, we use Eq. (j29|) and apply the same 
projectors, finding 

q 5 = nTr (v 5 T - ^-V 5 p) , (34) 



3n 



for the heat flux and 



P«*P> = - n TrV <a U p> . (35) 
5 

for the pressure deviator. With these expressions, the transport coefficients take the 
form: 

3c 2 k B 3 
kaw = — o — nr ' Vaw = -nk B Tr, haw = . (36) 
8 5 

Here, as in the case of the Marie model, we have restored the physical units. 
Note that the main difference between the transport coefficients, apart from different 
numerical prefactors, is that the ones calculated with the Anderson and Witting collision 
operator have a different dependence on the temperature than the ones calculated with 
the Marie model (since £ also depends on T). 

Considering a non-degenerate gas of relativistic particles in the ultrarelativistic 
regime, the number of particles density is given by n = k 2 B T 2 /2nc 2 h 2 . By taking this 
into account, we see from Fig. CD that the thermal conductivity km decreases with £ 3 , 
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Figure 1. Thermal conductivity k as a function of £. In this calculation we have 
set t = tm- The Marie coefficient is systematically higher than the Anderson- Witting 
one, and the ratio of the two grows at increasing l/£, i.e. at increasing temperature. 

while kaw with £ 2 . On the other hand, in Fig. [21 we can observe that the shear 
viscosity, displays the same qualitative behavior, t)m decreases with £ 4 while t\aw with 
£ 3 . In general, given any relativistic system, one can test which single relaxation time 
approximation, Marie or Anderson- Witting, better reproduces its behavior. We have 
considered only values l/£ > 1, since our calculations are valid only in this regime. 

An interesting calculation is to apply the corrections proposed by Takamoto pQ to 
the Marie model to account properly for the ultrarelativistic regime, £ — > 0, of a gas of 
particles. Although this work was developed in 3 + 1 dimensional space-time, we will 
follow a similar procedure for the case of 2 + 1 dimensions. To this purpose, we replace 
the relaxation time Tm by its average in momentum space, namely: 

mc f d?p 

TM = / —Frf^Trel = ^T rel , (37) 

n J p 

where r re i is now the effective relaxation time of the system and tm a simple parameter 
in the relativistic Boltzmann equation. By replacing this relation in the equations for 
the transport coefficients in the case of the Marie model, we obtain 

K re i = 3c 2 k B nT re l, 7] re i = 3nk B TT re i, /J, re l = . (38) 

Note that these transport coefficients carry the same dependence on temperature as in 
the case of the Anderson- Witting model. The numerical coefficients, though, are not 
the same. 

5. Conclusions and Discussions 

In this work, we have calculated the transport coefficients, namely the bulk and shear 
viscosities, and the thermal conductivity of a two dimensional ultra-relativistic system, 
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Figure 2. Shear viscosity 77 as a function of £. In this calculation we have set 
t = tm- Like for the case of Figure 1, the Marie coefficient is systematically higher 
than the Anderson- Witting one, and the ratio of the two grows at increasing l/£, i.e. 
at increasing temperature. 



using two different forms of the collision operator. The first one is based on the 
Marie model and the second one on the Anderson Witting approach. Depending on 
the approach, we have to satisfy the Eckart or the Landau-Lifshitz decompositions, 
respectively. This leads to different expressions for the transport equations and third 
order moment of the distribution function. 

We have found that the bulk viscosity of the ultrarelativistic system disappears as 
a consequence of the choice of the two-dimensional ultra-relativistic equation of state, 
which imposes a constraint on the trace of the momentum-energy tensor. This is the 
same behavior observed for the three dimensional case. By analyzing the transport 
coefficients for the case of an ultrarelativistic gas of particles, we have found that the 
thermal conductivity decreases with £ 3 and £ 2 , for the case of Marie and Anderson- 
Witting, respectively. The shear viscosity presents the same qualitative behavior, 
decreasing with £ 4 and £ 3 for both models, respectively. Therefore, the Marie model 
transport coefficients always decrease faster than the ones based on the Anderson- 
Witting model. In a more general relativistic system, by knowing this difference, one 
could select which one of the two is better suited to describe its dynamics evolution. 
In addition, following the work by Takamoto [1] , we have modified the two-dimensional 
transport coefficients for the case of the Marie model, in such a way as to make it 
suitable for a gas of ultra-relativistic particles. With such modification, the functional 
dependence of the transport coefficients on the temperature becomes the same as for 
the Anderson- Witting model, although with different numerical coefficients. 

It is known that transport coefficients in 2d are formally infrared divergent, hence 
their size and gradient dependence must be taken with some caution in practical 
applications [T8l [19]. The investigation of these issues in the relativistic framework 



Ultrarelativistic Transport Coefficients in Two Dimensions 



12 



makes an interesting object of future research. 

The results presented in this paper can be applied to a variety of ultrarelativistic 
systems, e.g. graphene, plasma jets and others. The method is not limited to ultra- 
relativistic gases of particles, and it extends to any statistical system obeying relativistic 
Boltzmann-like equations. 

Acknowledgments 

We acknowledge financial support from the European Research Council (ERC) 
Advanced Grant 319968-FlowCCS. I.V.K. was supported by the ERC Grant ELBM. 

Appendix A. Moments of the Equilibrium Distribution 

The moments of the equilibrium distribution for a two-dimensional ultrarelativistic 
system that satisfies the Maxwell Juttner distribution are given by, 

Ng = nU a , (A.l) 
Tf = -nTrf p + ZnTU a U p , (A.2) 

Tjf 7 = -3nT 2 {r] a/3 U^ + + rfrU a ) 

+ lbnT 2 U a U p W , (A.3) 

T^ 5 = 3nT 3 (rj a % lS + rj a ^ s + rf p rf & ) 

- IhnT^^U^U 5 + n^U^U 6 + rf p U a U s 
+ r) aS WUP + r]^U^U a + r] S ^U a W) 

+ l05nT 3 U a UPlT f U 5 , (A.4) 

+ U a {r} eP rf 5 + rfV 5 + TfV 4 ) 
+ U p (ri ae rf & + r] a ^r] eS + T? € T) aS ) 
+ U p {ji M rp s + r^V 5 + V le V aS ) 

+ ^{rj^rf 6 + 77 a V £ + ^V')] 

- IQhnT^^W^U' + V a ^U p U 5 U e 
+ r, a5 U^WU e + r,^U a U^U e 

+ rfPU a WU e + V ^U a U 5 U e 

+ ri ae U^WU 5 + r]^U a WU 5 + r]^U a U^U 5 ) 

+ MbnT 4 U a U p WU 5 U e . (A.5) 
To obtain the moments, we have considered the ultrarelativistic regime, £ <C 1. 
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